Effect of Impurities on the Superheating field of Type II superconductors 
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We consider the effect of nonmagnetic and magnetic impurities on the superheating field Hs in a 
type-II superconductor. We solved the Eilenberger equations, which take into account the nonlinear 
pairbreaking of Meissner screening currents, and calculated Hs{T) for arbitrary temperatures and 
impurity concentrations in a single-band s-wave superconductor with a large Ginzburg-Landau pa- 
rameter. At low temperatures nonmagnetic impurities suppress a weak maximum in Hs{T) which 
has been predicted for the clean limit, resulting instead in a maximum of Hs as a function of impurity 
concentration in a moderately clean limit. It is shown that nonmagnetic impurities weakly affect 
even in the dirty limit, while magnetic impurities suppress both Hs and the critical temperature Tc- 
The density of quasiparticles states A'^(e) is strongly affected by an interplay of impurity scattering 
and current pairbreaking. We show that a clean superconductor at H — Hs is in a gapless state, but 
a quasiparticle gap eg in N{e) at H = Hs appears as the concentration of nonmagnetic impurities 
increases. As the nonmagnetic scattering rate a increases above Oc = 0.36, the quasiparticle gap 
eg (a) at H = Hs increases, approaching eg ~ 0.32Ao in the dirty limit a 3> 1, where Ao is the 
superconducting gap parameter at zero field. The effects of impurities on Hs can be essential for 
the nonlinear surface resistance and superconductivity breakdown by strong RF fields. 



I. INTRODUCTION 



Type -II superconductors are in the Meissner state if 
the applied magnetic field H is smaller than the lower 
critical magnetic field Hd above which the Gibbs free 
energy of a vortex becomes negative. However, the Meiss- 
ner state can remain metastable at higher magnetic fields, 
H > Hci up to the superheating field Hg at which the 
Bean-Livingston surface barrier^ for penetration of vor- 
tices disappears and the Meissner screening currents at 
the surface become unstable with respect to small pertur- 
bations of the order parameter. The field Hg is therefore 
the maximum magnetic field at which a type-II super- 
conductor can remain in a true non-dissipative state not 
altered by dissipative motion of vortices. In addition to 
the fundamental interest in what limits Hs in supercon- 
ductors, there is a strong interest in the physics of Hs 
promoted by recent advances in superconducting Nb cav- 
ities for particle accelerators in which the quality factors 
Q - 10^° - 10" at 2K and 2GHz were observed up to 
the surface RF fields close to the thermodynamic critical 
field He ~ 200 mT of Nb,^ Here the RF frequency is well 
below the gap frequency Ag/h ~ 400 GHz for Nb where 
Aq is the superconducting gap, so the current pairbreak- 
ing in the cavities at low temperatures may occur under 
quasi-static conditions. 

Calculations of the superheating field Hg and the de- 
pairing current density Jc have a long history starting 
from the pioneering works by Ginzburg'^, de-Gennes4 
and Matricon and Saint-James"' who used the Ginzburg- 
Landau (GL) theory. It was shown that as the applied 
field H reaches Hs, the Meissner screening current den- 
sity at the surface becomes of the order of Jc ~ cHc/4:7tX, 
which makes the superconducting state unstable. The 
value of Hs in the limits of large and small GL parame- 



ter K is given hy^ 
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Hs = {V5/S)Hc, 



Hs = 2 



- 9-1/4^-1/2 



He, 



K > 1, 
K < 1. 



(1) 

(2) 



Here the factor -\/5/3 in Eq. ([T]) reflects the suppression of 
the order parameter by the Meissner currents in the local 
limit K ^ 1, while the enhancement of Hg by the factor 
fQj- K ^ 1 results from the proximity effect reduc- 
tion of the Meissner pairbreaking localized in a narrow 
surface layer of thickness A <C where A is the London 
penetration depth and ^ is the coherence length. The ex- 
tensive calculations of Hs based on the GL equations^"— 
have shown that, for k ^ 1, the Meissner state becomes 
absolutely unstable with respect to small 2D perturba- 
tions of current and order parameter with the wavelength 
~ (C"^A)^/^ along the surface and decaying over the length 
~ \/Aj perpendicular to the surface. Such perturbations 
describe the initial stage of penetration of vortex rows 
with the period ^ {(j)o/HsY^'^ corresponding to the equi- 
librium vortex lattice at the field H = Hg. 

Unlike Hs{T) in the GL region near Tc, the behavior of 
Hs{T) at low temperatures is not well understood, not 
least because the calculation of Hs{T) requires solving 
the nonlinear Gor'kov or Eilenberger equationsi^ which 
take into account pairbreaking effect of Meissner cur- 
rents. Manifestations of pairbreaking effects in a clean 
superconductor at low temperatures can differ from the 
GL results, as was shown long ago by Parmenter— and 
Bardeeni^. The first calculation of Hs (T) for the entire 
temperature range < T < Tc in the clean limit and k — 
oo was done by Galaikc^^ who obtained Hs — O.SAHc at 
T ^ 0, and Hs = {V5/S)Ha = 0.745i?c at T ^ Tc. 
Catelani and Sethna^^ solved the Eilenberger equations 
to calculate the temperature dependence of Hs{T) for 
< T < Tc in the clean limit for /t — > cx) and found 
a maximum in Hs{T) at low T. Such a non-monotonic 
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temperature dependence of Hg (T) shows that the behav- 
ior of Hg(T) at low T can hardly be extrapolated from 
the GL results near Tc- 

The effect of impurities on Hg{T) outside the GL re- 
gion has not been addressed. This problem is of inter- 
est because the clean limit in s-wave superconductors at 
T <C Tc is a rather singular case: for T = 0, the Meissner 
currents do not affect the superfluid density ng until the 
superfluid velocity Vg — J/ugC reaches the critical value, 
Vs — Vc — Aqq/pp where is the Fermi momentum 
and Ago — Aq (0) is the modulus of the order parameter 
vl^ = A exp{iip) at zero superfluid velocity and T = O i^'^i^^ 
For Vg > Vc, the gap in the quasiparticle spectrum disap- 
pears and ng(vg) rapidly drops to zero in a narrow region 
Vc < Vg < l.OSwc'^'^ Thus, unlike the d-wave supercon- 
ductors, the s-wave superconductors at T ^ do not 
exhibit the nonlinear Meissner effect caused by the de- 
pendence of ng{J) on the current densityi^i"— , and the 
superheating field is reached at the superfluid velocity 
Vg > Vc, which corresponds to a gapless state^^. The 
latter has important consequences for the low-frequency 
(w <C Aq) impedance of clean superconductors because 
aX H = Hg = O.SAHc the surface resistance Rg becomes 
of the order of Rg in the normal state, unlike the expo- 
nentially small Rg ~ exp(— Ao/fcsT) in a fully gapped 
state at 7J < HgM'^ The effect of the RF Meissner cur- 
rents on the quasiparticle spectrum can result in a strong 
dependence of the surface resistance on the RF field am- 
plitude with H > THc/Aq.^^- 

For k:s> 1, the field Hg at which the gap in the quasi- 
particle density of states N{e) closes at T = can be cal- 
culated from the London equation Hg = AttXJc/c where 
Jc = nevc, A = (mc^ /Anne^y^^ , n is the total electron 
density, — e is the electron charge, m is the band effec- 
tive mass, and c is the speed of light. Here the linear 
London screening of H(x) = HQe~^^^ remains valid as 
long as the superfluid density is independent of J, that 
is, the Meissner current density J = {c/4TT)dH/dx is 
smaller than Jc- This yields Hg = cAqq/cXvp, which 

1/2 

can be expressed in terms of He — (47riV(0)) Aqo, 
where iV(0) = w?vf /'^T^'^fi^ is the density of states per 
one spin orientation for an isotropic parabolic band, and 
vf — {3TT^ny/^h/m. Hence, 

Hg^{2/sy/^Hc~0M6Hc. (3) 

Since Hg — O.SAHc,^- the gapless state in the clean, large- 
K limit occurs in a narrow field range, 0.97iJs < H < Hg. 

The goal of this work is to address the pairbreaking 
effect of impurities on Hg. Indeed, while nonmagnetic 
impurities do not affect Tc, Aq and He at zero super- 
fluid velocity2Zi2S, these impurities become pairbreakers 
in the current-carrying state^^, which, for example, man- 
ifests itself in the nonlinear Meissner effecti^i"— . Thus, 
the extent to which nonmagnetic impurities would affect 
Hg needs to be understood. The effects of impurities 
on the quasiparticle density of state N{e) at H = Hg 
and the conditions under which impurities can restore 



the gap in N{e) a.t H = Hg are important for the un- 
derstanding of the nonlinear surface resistance at high 
fields and the limits of superconductivity breakdown un- 
der low frequency RF fields. We will also consider the 
superheating field in thin film multilayers consisting of 
alternating superconducting and dielectric layers thinner 
then A. The parallel Hd of such multilayers is greatly 
enhanced, which enables probing the pairbreaking limits 
in parallel magnetic fields. It was suggested to use these 
multilayer coatings to increase the RF breakdown fields 
of superconducting cavities in particle accelerators The 
paper is organized as follows. 

In Section II we solve the Eilenberger equations for a 
superconductor with uniform current and analyzed the 
nonlinear dependencies of A and J on the superfiuid ve- 
locity in the presence of impurities. In Section HI we 
obtained equations for Hg{T) in the entire temperature 
range < T < Tc for k ^ 1 and arbitrary concentrations 
of nonmagnetic and magnetic impurities. We solved the 
equations for Hg numerically and analyzed the depen- 
dencies of Hg on the magnetic and nonmagnetic scat- 
tering rates. The range of parameters in which Hg can 
be optimized by varying the impurity concentration was 
found. In Section IV we consider the effect of impurities 
on the quasiparticle density of states a.t < H < Hg, 
particularly the emergence of the gap in N{e) as the con- 
centration of nonmagnetic impurities increases. Section 
V is devoted to the calculation of nonlinear screening of 
strong DC field and Hg in multilayers. In Section VI we 
discuss our results and their implications for a nonlinear 
surface resistance at high RF fields and superconductiv- 
ity breakdown in cavities. 



II. THEORY 
A. Eilenberger Equations 

In this work we use the Eilenberger equations for a 
single-band s-wave superconductor with magnetic and 
nonmagnetic impuritiesJ^ 

(2w„ + v2?)/ = 2M'g+— g(/>-— /(.g), (4) 
r_ r+ 

where ^' is the order parameter, v is the Fermi velocity, 
= 7rT(2n -|- 1) are Matsubara frequencies, T is the 
temperature, I? = V + 27riA/(/)o, is the flux quantum, 
and A is the vector potential. We use the units for which 
h = ks — 1 unless stated otherwise. The time constants 
T+ and T_ are defined by 



where r is the electron scattering time on nonmagnetic 
impurities, and is the spin-fiip scattering time on mag- 
netic impurities. The angular brackets (• • • ) = Jg^ d^kp 
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mean angular averaging over the Fermi surface. The qua- 
siclassical Green functions /(v,r,w„) and (7(v, r,cj„) are 
normalized by 

.9' + //'=!• (6) 

where /'''(v,r, w„) = /*(— v,r,a;„), and the asterisk 
means complex conjugation. The self-consistency equa- 
tion for the superconducting order parameter ^(r) = 
A exp(iip) is given by 

* = 27rriV(0)|y| ^ (/), (7) 

where A'^(O) is the normal density of states per one spin 
at the Fermi surface, V is the BCS coupling constant 
and the cut-off frequency uju is of the order of Debye 
frequency. Here V can be expressed in terms of the crit- 
ical temperature Tc of a superconductor without non- 
magnetic impurities: \V\~^ — N{0)\n{2ujuj/TrTc) where 
7 = 1.78. Eqs. (H])-® are supplemented by the Maxwell 
equation for A and the supercurrent density J: 

An 

V X V X A = —J (8) 

c 

00 

J = -47rreiV(0) Im ^ (jrv). (9) 

LJ„=0 

B. Geometry and assumptions 



the coordinates. Equations (ITOl) and ([TT|) approximate 
the solution necessary for the calculation of Hs, taking 
into account the difference in characteristic lengthscales 
in an extreme type-II superconductor with k 3> 1. In this 
case Hs is reached as the current density at the surface 
becomes of the order of the depairing current density, 
which corresponds to the wavelength 27r/g ~ Slow 
decrease of the Meissner screening currents over the Lon- 
don penetration depth increases Hs which is now limited 
by perturbations S3{x,y) and 5f{uJn,Q,x,y) oscillating 
on the scale 27r/fc ~ {Sj^^Y^'^ along the x— axis and de- 
caying over ^ ii^Y^"^ along the y— axis2r— lii. Thus, Eqs. 
pO| and (ITT]) give an asymptotically exact solution for 
the lower bound of Hs in the limit of k — >■ 00 which will be 
addressed in this paper. The effect of London screening 
resulting in the finite-fc instability produces small correc- 
tions in Hs = i?c(\/5/3)(l + l/%/2^) in the GL region^iii. 

Equations pUj) and ([TT|) can also be applicable for a 
wider range of k in a multilayer system in which super- 
conducting layers are separated by thin dielectric layers 
as shown in Fig. [T] For thin superconducting layers of 
thickness ci <C A, the Meissner current is nearly uniform 
across the film, while for d < \/X^, the vortex instability 
is suppressed as the perpendicular components of J{x, y) 
cannot cross the dielectric layers.— In this case a uniform 
pairbreaking instability at H = Hs develops, for which 
equations PH)) and are exact solutions of the Eilen- 
berger equations. These solutions will be used to calcu- 
late nonlinear screening of magnetic field in multilayers. 



To study the stability of the Meissner state with re- 
spect to small perturbations, we consider a planar type-II 
superconductor occupying the region a: > with mag- 
netic field H applied along the z-axis. For a large-K 
superconductor, the previous calculations based on the 
Glr^^^^ and Eilenberger— theories have shown that the 
instability at H = Hg is driven by coupled fluctuations of 
(5^(x, y) and 53{x,y) that rapidly oscillate on the scale 
~ ^K^/** parallel to the surface, and decay on a longer 
length ~ ^K^/^ perpendicular to the surface. Taking 
these inhomogeneous unstable modes into account is es- 
sential for the calculation of small corrections ^ H^k^^^^ 
to Hs- However, the main contribution to Hg He is 
determined by the condition that the superfluid velocity 
v{Q) at the surface reaches the critical pairbreaking value 
Vc for which J(v) is maximum and the superconducting 
state becomes unstable. It is the same condition that 
defines the uniform depairing current density^™. 

For uniform current flow, the solution of the Eilen- 
berger equation can be sought in the form: 



*(r) = Ae'«^ 



(10) 
(11) 



where q is the wave vector of the condensate, where 9 
is the angle between the Fermi velocity and the current, 
and the amplitudes A and f{0,uj„) are independent of 



o 



t 







FIG. 1: Schematic diagram of a multi-layered superconduc- 
tor: Superconducting layers (gray) separated by insulating 
layers (dark lines). Applied magnetic field applied along the 
z-direction induces currents fiowing along j/-direction. 
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C. Solution of the Eilenberger equations 

Substituting Eqs. ^ and (dH into Eq. (g]) gives the 
following algebraic equations for the amplitudes g[6,ujn) 
and f{0,ujn): 



9 1 ^ (.9) ^ . 

7 ^ 2A +(/)/._ 1^-" + -+—' 



(12) 



where we assumed a spherical Fermi surface for which 
(/) = i f{6) sm9d9, 9 is the angle between the Fermi 
velocity and the direction of current. The gauge-invariant 
wave vector u = q + 2'!TAy/(j)o is proportional to the su- 
perfluid velocity Vg — hu/2m, where m is the effective 
electron mass. Solving Eqs. (|6]) and (fT2|) yields 



/ = 



^^52 + {a + icos9f 



a + I cos ( 



y^52 + (a + icos6l)^ 
where a and b depend on (g) and (/): 

a ^ (2w„ + (g)/T+)/uv, 
b = {2A+ if )/t^)/uv. 



(13) 



(14) 



(15) 
(16) 



Integrating Eqs. ([T3l) and (fT4|) over 6*, we obtain two self- 
consistent equations for (/) = —b Imsinh^^ ((a + i)/b) 

and (g) = Im^b'^ + {a + i)"^. The equations for (/) and 
((?) can be recasted in the following form: 



(/>-&tan-iM 
a 



(17) 

{g)^ + {gf{a^ + b^-l)-a^ = 0. (18) 
It is convenient to introduce a new variable X such that 



X = 



(g) 



2Ar_ + (/) ' 

2uJnT+ 

UVT+ cot X ~ 1 



(19) 
(20) 



Substituting this into Eqs. (ITTl) and gives the fol- 
lowing equation for X 



A sin Xt^ 



uvT- — X J \ uvtj^ cot X — \ 



1 
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(21) 



The gap equation and the current density J = |J| in Eqs. 
([7]) and ([nj can be expressed in terms of X as follows: 



2Xt_ 



= 0, (22) 



J = 4^re«7V(0)r!A2f;^^-^^. (23) 



^ (UWT_ - X) 



Equations ([T9])-(l23l) imphcitly define (/), (5), and A 
as functions of J in a superconductor with arbitrary con- 
centration of nonmagnetic and magnetic impurities in the 
entire range < T < T^. Generally, Eqs. (IT1|)-(123]) can 
be solved only numerically, but analytical solutions can 
be obtained in some limiting cases. 

To make the effect of impurities on Hg more clear, 
we first discuss the instructive case of a clean supercon- 
ductor (T4.,r_) — ooJ^ii^ii^ In the absence of current 
u ^ 0, Eqs. (HH) and ^ yield the Gor'kov result, 
(/) = A/VA2 +w2. Solving Eqs. (HH) and ^ for a 
current-carrying state, u > 0, gives: 



(24) 

^^ 1/2 

(25) 



4^2 



A2 



(/> = 


2A 1 r- 


tan 


uv 


-1 + 


( /u^v^ A2 


I V 4^2 UJ^ 



The dependence of A on u and T is then determined by 
Eqs. 0, (IMl) and 

Let us discuss the clean limit at T = in more detail. 
For T — 0, the summation in Eqs. ([^ and (^5]) can be 
replaced by integration. Substituting Eq. (1^^ in those 
equations gives 



.18.30 



In- 



A, 



00 



0, 



u; < 1 



cosh W + y/l — l/iy2, W > 1 



(26) 



where w = uw/2A, A is function of u, and Aqo is the 
order parameter at u = and T — 0. Calculating J(u) 
in Eq. ([23]) yields 



Jo -\l-(l-l/z.2)3/2, 



W < 1 
W > 1. 



(27) 



Here the current density 



Jo = —eN{0)v^u = envs- 
3 

corresponds to superconducting flow of all electrons. 

The calculated dependencies of J{u) and A(u) on u 
for a clean superconductor at T — ^ are shown in Fig. [2] 
The current J = envg increases linearly as m = 2mvs/h 
increases up to u = 2Aoo/w, then J{u) becomes nonlinear 
and reaches the maximum value Jc at the pair-breaking 
momentum Us — 2.059Aoo/w, and then drops to zero. At 
the same time, the gap A(u) = Aoo remains unaffected 
by current at < u < 2Aoo/w, but rapidly drops to zero 
for u > 2 Aqo /v. Here Ag at u = Us is slightly smaller 
then Aoo. 

A clean s-wave superconductor at T = exhibits no 
current pairbreaking in the broad range of < u < 
2Aoo/i' where the linear London electrodynamics is ap- 
plicable. For finite temperatures, this anomalous fea- 
ture no longer holds, as illustrated by Fig. 2 which 
shows J{u,T) and A(u,T) calculated numerically from 
Eqs. ©-(IMl) for T = 0.5Tc. Here A(u) decreases 



5 



1.2 
1 

0.8 

1^ 0.6 
0.4 
0.2 




1 


1 1 1 


(a) 




It c-. 


\t=o 




t=0.5'X^ 




1 


1 1 1 





0.5 



1 1.5 

n-u/Aoo 



2.5 



1.2 
1 

0.8 
0.6 
0.4 
0.2 




1 


1 1 1 








\ t=0 




t=0.5'''\ 






1 1 1 





0.5 



1 1.5 



2.5 



FIG. 2: (a) Current density and (6) gap as functions of u 
at T -> (solid line) and at T = 0.5Tc (dashed line). For 
T = 0, 3(u) increases linearly with w, reaches Jc at Us and 
decreases to zero. The gap stays constant for u < 2Aoo/f 
and then drops. The arrow shows the location of the depar- 
ing momentum its ~ 2.059Aoo/i' where current reaches the 
critical current Jc- At T = 0.5Tc, Us is shifted to lower u. 



as u increases, and J{u) deviates from linearity for u 
smaller than Aqq/v, similar to the nonlinearity of J(u) cx 
m(1 — u^^^) due to current pairbreaking in the GL theory. 
Such a difference in the behaviors of A(w) at T = and 
T ~ Tc is due to thermally-activated quasiparticles that 
are negligible at T <C Tc but become essential at higher 
temperatures. In any case, J{u) reaches the maximum 
at the critical superfluid momentum u ~ v/Aqq ^ 
where = Sw/ttAoo is the superconducting coherence 
length in the clean limit. The maximum of J(u) at 
the critical momentum Us defines the superheating field 
which will be calculated in the next section. 

Impurities do not change qualitatively the behavior of 
J{u) shown in Fig. |31 but they increase the critical mo- 
mentum Us. We calculated the effect of nonmagnetic 
impurities on at T — >■ by solving the full Eqs. (PT|). 

and The results are shown in Fig. [3] which 

displays function of the scattering parameter 

a = I/Aoqt — TT^o/^ where £ = vt is the mean free 
path. The critical momentum Us increases as a increases, 
approaching the square root dependence Us oc ^/a for 
a > 1, consistent with the general result, Ug ^ 1 /£,, where 
^ ~ {^o^y^^ in the dirty limit. The dependence of Us{a) 
can be described by the interpolation formula 



(l.053(a-|-0.655) 



1/2 



which approximates the calculated u, 
to an accuracy better than 1%. 



1.146)Aoo/i^, (28) 
a) for < a < 20 



Numerical solutions of the equation for A at T = 
show that in the presence of impurities A(m) decreases 
as u increases even if ^ 1. As a result, nonmag- 
netic impurities become pairbreakers, reducing the su- 
perfluid density as J increases,^- This manifests itself in 
a nonlinear Meissner effect even at T = 0.^^ Shown in 
Fig. [3] are the calculations of Ag at the depairing mo- 
mentum Us as a function of the scattering parameter a 
a,t T — 0. Here As decreases as a increases, approaching 
As (a — > oo) = 0.79Aoo in the dirty limit. 



CO 

;3 




FIG. 3: (a) The pairbreaking momentum Us{a) as a function 
of the nonmagnetic scattering rate a at T = 0. (b) at 
the critical momentum u = Us as a function of a. Here As 
approaches 0.79Aoo for a ^ 1. 



III. SUPERHEATING FIELD 

We calculate Hs for a superconductor occupying the 
half-space a; > by solving the Maxwell equation 



dB At: 
dx 



-J, 



(29) 



where the field is applied along the z— axis, so both the 
current density J{u) and the vector potential Ay have 
only the y— components. Here J{u) — —c5J-/5Ay in Eq. 
([9|) can be obtained by varying the free energy densityi^ 



T °° r A^ 

F = A2iV(0) In — + 27rAr(0)r V <^ 2A(/)- 

n— ^ 

2wn{{g) - 1) - iuv{gcose) - M ^■'^ ' 



(30) 
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where / and g are the sohition of Eqs. (|4]) and ([6]). For 
K ^ 1, the amphtudes A and / vary slowly on the scale 
of Us ~ 1/^ so the gradient terms V/ in are negligi- 
ble. In this case the variational derivative J — —c5F/5A 
becomes the partial derivative, J = —cdjFjdA. More- 
over, since dTjdA = dTjdA + [pT ldl\){dt\ldA) + 
{dTldj){dI^ldj) where dTjdlA = dTjd] = in equi- 
librium, we have J = —cdTjdA. 

Multiplying both sides of Eq. (gHl) \iy B ^ dAj dx and 
integrating from x = to a; = cx), we obtain 



' dS^ ^ [°^dTdA^ 
B^da; = -47r / — ^— da;. 



dx 



dA dx 



(31) 



Here the boundary conditions are: i? and T ^• 
-Hll^-K at X — >■ oo, and B ^ H x — Thus, Eq. 
([5T|) reduces to: 



ij2 = ij2 + 87rJ"[M(0)], 



(32) 



where i?c is the thermodynamic critical field and J^[u(0)] 
and u(0) are the free energy density and the momentum 
of condensate at the surface, respectively. Eq. (15^ thus 
determines the dependence of u(0) on the applied field. 

As the magnetic field B increases, J(0) and u(0) at 
the surface increase. However, as u(0) reaches Ms at 
which J(u) is maximum, further increase of B does not 
cause any increase of screening current density, making 
the Meissner state absolutely unstable with respect to 
small perturbations of the order parameter and initiat- 
ing penetration of vortices. The condition u(0) = Ug in 
Eq. ([5^ defines the superheating field: 



Bl{T) = Bl{T) + %^T{u,,T). 



(33) 



Substituting the gap equation ([7]) and 2i{gcosd) ~ 
a{{g) + i/(.g)) - b{f) to Eq. (1201) yields J- = 
27rA^(0)TX;^„>oW„(2 - (g) - l/{g)). Hence, we obtain 
the final expression for Bs'. 



B^ 



Bl + IQtt'^N{Q)T 



n=0 



{9s 



(34) 



Here the index s means that the function {g) is calculated 
at the critical momentum u = u^. The thermodynamic 
critical field Be is given by 



B^ = 16tt'^N{Q)tY^ 



- 2un } (35) 



where Ao(r) is the superconducting order parameter at 
u = 0. In the clean limit, Eq. (p4|) reduces to: 



B^ = 16n'^TN{0) 



(36) 



where z{u) is given by Eq. (|24p 



Equations ([341) and ([35]) combined with Eqs. ([2T]) - ([23)) 
define Bg as a function of temperature, and concentra- 
tions of nonmagnetic and magnetic impurities. Calcu- 
lations of Bs involves solving coupled Eqs. (HU, (|22p . 
p3| to obtain A and J as functions of u and then find- 
ing self-consistently the depairting momentum Us from 
the condition dJ/du = 0. The results of these calcula- 
tions of Bs for different concentrations of impurities and 
temperatures are presented below. 



A. Effect of nonmagnetic impurities 

We quantify the effect of nonmagnetic impurities on 
Bs by the dimensionless scattering rate 



(37) 
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where = ^i'/'^'Aqo is the clean limit coherence 
length, and ^ = ut is the mean free path. For u = 
0, impurities do not affect A, Be and Tc of s-wave 
superconductors22i^j but in a current state both A and 
Us are affected by impuritiesi^, as shown in Fig. [3] 

We first consider the effect of nonmagnetic impuri- 
ties on the temperature dependence of Bs{T). Shown 
in Fig. m are the superheating fields calculated for the 
clean (a = 0), moderately dirty (a = 1) and dirty 
(a = 10) limit. The first conclusion apparent from Fig. 2] 
is that nonmagnetic impurities have a rather weak effect 
on Bs, despite their pairbreaking nature in the presence 
of current. Near Tc the BsiT) curves coincide, approach- 
ing Bs{T) = {VE/3)Bc predicted by the GL theory, 
where Be is independent of a according to the Ander- 
son theorent^. At lower temperatures impurities result 
in different behaviors of Bs{T) for different a. For in- 
stance, one of the features of the clean limit at low tem- 
peratures is a maximum in Bs{T) at T = 0.04Tc pointed 
out by Catelani and Sethna.— As follows from Fig. |4j 
this maximum disappears in the moderately dirty limit 
{a = I), as well as the dirty limit (a = 10). As will 
be shown below, this change in the behaviors of Bs{T) 
at low temperatures results from the impurity-induced 
change in the quasiparticle density of states at B — Bs 
from the gapless state in the clean limit to a gapped state 
for a > ac- 

We now turn to the dependence of Bs on the nonmag- 
netic scattering rate a. Shown in Fig. [S]are Bs{a) cal- 
culated for different temperatures. These Bs{a) curves 
exhibit surprising maxima at low temperatures, resulting 
in a ~ 4.2% enhancement of Bs{a) at T = as compared 
to the clean limit. Here the arrows in the figure mark 
the maxima in Bs, the location of which depends non- 
monotonically on temperature as described in the cap- 
tion to Fig.O Thus, there is an optimum mean free path 
^maa = 5.3^0 at T = for which the superheating field 
is maximum. At higher temperatures, T > 0.3rc, the 
maximum in Bs{a) disappears. 

The results of this section show that the effect of non- 
magnetic impurities on Bs is most pronounced at low 



7 



0.8 




0.9 

^ 0.8 



(a) 




a = 
1 
10 



0.1 0.2 



0.2 0.4 0.6 

T/T, 



0.8 



1 1 


1 1 

(b) 


1 1 


1 1 ^ 



0.2 



0.4 



0.6 



0.8 



FIG. 4: Temperature dependencies of Hs calculated for a = 
0,1 and 10. (a) H^T) in the units of Hc{0) and (b) the 
ratio Hs(T)/Hc{T) where Hc{T) is calculated from Eq.(|35}. 
At low T, nonmagnetic impurities result in nonmonotonic de- 
pendence of Hs on a, eliminating the maximum in Hs{T) at 
T = 0.04Tc for the clean limit, as shown in the inset. The 
difference between Hs{T) diminishes at higher T where the 
Hs{T) curves approach the GL result given by Eq. ([l}. 



temperatures where impurities can eliminate the maxi- 
mum in the temperature dependence of Hg characteristic 
of the clean limit. At the same time, impurities can cause 
a new maximum in Hs as a function of the scattering rate 
a at low temperatures, although the overall effect of non- 
magnetic scattering on Hg turns out to be comparatively 
weak. This results from two opposite effects that nearly 
cancel each other: the increase of Hg due to the increase 
of the pairbreaking momentum Us shown in Fig. [3] and 
the decrease of superfluid density as a increases. Here 
Hg roughly scales like the thermodynamic critical field 
He in which these opposite trends cancel out exactly, as 
prescribed by the Anderson theorem. Yet our results 
show that this cancelation is not exact for Hg which de- 
pends weakly on the scattering rate a, as illustrated by 
Fig. [B] which shows the calculated ffs(0, a) at T = 0. 
For instance Hg « 0.82Hc at a = 10 is rather close to 
Hs = O.SAHn in the clean limit. 




FIG. 5: Enhancement of Hs by nonmagnetic impurities at 
low temperature where Hs{a,T) is normalized to Hso{T) for 
a clean superconductor. Here the position of the maximum 
OmaxiT) in Hs{a) depends on temperature: a,nax ~ 0.59 for 
T = 0, Umax = 0.91 for T = ClTc, and amax = 0.53 for T = 
0.2Tc. At higher temperatures, the maximum disappears. 





0.88 r 




0.87 ■ 




0.86 






0.85 - 




0.84 - 




0.83 - 




0.82 - 




0.81 - 




FIG. 6: Dependence of Hs at T — on the nonmagnetic 
scattering rate a. 



B. Effect of magnetic impurities 

Unlike nonmagnetic impurities, magnetic impurities 
strongly suppress Tc,^* resulting in a significant reduc- 
tion of Hs as well. We quantify this effect by the dimen- 
sionless magnetic scattering rate am similar to that was 
used above for nonmagnetic impurities: 



(38) 



where and im = are the spin-flip scattering time 
and mean free path, respectively. 

Shown in Fig. [7] are the temperature dependencies 
Hs{T) calculated for a = and different values of am- 
The superheating field is significantly reduced by mag- 
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netic scattering, the Hs suppression is stronger at higher 
temperatures. One can see that for — 0.1, which cor- 
responds to £m ^ 30^0: the superheating field is roughly 
reduced by half as compared to the case of am ~ 0. 




T/T,(0) 

FIG. 7: Temperature dependencies of Hs{T) calculated for 
a = and different magnetic scattering rates a,n- Here Tc(0) 
is the critical temperature for clean superconductors. 



Similar to the case of am = 0, the effect of nonmagnetic 
impurities at am > is most noticeable at low temper- 
atures, where it remains rather weak even for very dirty 
superconductors with a ~ 100. 




FIG. 9: Hs{am) calculated for T = 0.2rc(0) and different 
nonmagnetic scattering rates a. 



Figure |S] shows the reduction Hg{am) by magnetic im- 
purities at a given temperature. Unlike the effect of 
nonmagnetic impurities which is most pronounced at 
low temperatures, the effect of magnetic impurities is 
stronger at higher temperature. Here Hs{am) decreases 
as am decreases, vanishing at a critical scattering rate 

OLmc{T). 




FIG. 8: Dependencies of Hs on am calculated for at T = 
0.2rc(0) and T = 0.57^(0) Insert shows the dependence of Tc 
on oim calculated from the Abrikosov-Gor'kov theory^S 

The combined effect of magnetic and nonmagnetic im- 
purities on Hg is shown in Fig. |9l From the curves 
Hs{am) calculated for different a at T = 0.2Tc(0), we 
can see that there is practically no interplay between the 
effects of magnetic and nonmagnetic scattering on Hg. 



IV. QUASIPARTICLE DENSITY OF STATES IN 
THE MEISSNER STATE 



As was discussed in the Introduction, the effect of cur- 
rent on the quasiparticle density of states N{e) in the 
Meissner state is rather nontrivial because a clean su- 
perconductor a,t H = Hs is in a gapless state. The ap- 
pearance of quasiparticle at the Fermi level at Hs can 
have important consequences for the breakdown of su- 
perconductivity under strong RF fields, as will be dis- 
cussed below. This feature of the s-wave clean limit re- 
sults from the spectrum of quasi-particles in a uniform 
current-carrying stati 



,13.14. 



Ek = + v^(.P - PfY + ftu • v/2 



(39) 



where the last term describes the effect of current with 
the superfluid velocity Vs = hu/2m. The quasiparti- 
cle gap Eg corresponds to the minimum of Ek, giving 
eg — Aq for u — 0. However, in the presence of cur- 
rent, Ao(r) is no longer the gap in the quasiparticle 
spectrum. For instance, Aqo at T = remains inde- 
pendent of J in the entire interval < u < 2Aoo/fti', as 
shown in Fig. [21 while the anisotropic quasi-particle gap 
eg{u) — |Aoo| + ■^huv coaip, obtained by setting p = pf 
in Eq. (p9)) . depends on the angle (p between the cur- 
rent and quasiparticle momentum. The minimum gap 
Sg = I Aqo I — huv/2 vanishes at the condensate momen- 
tum u > Ug — 2Aoo/ftt' slightly smaller than the pair- 
breaking momentum Us corresponding to H = Hs- The 
question is then how the density of states at H = Hs is 
altered by impurity scattering. To address this issue we 
use the approach developed by Makii^ii^ and Fulde^^^ to 
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calculate N{e, u) as functions of energy e and the super- 
fluid velocity. For the sake of simplicity, we only consider 
here the effect of nonmagnetic impurities on A^(e, u). 

The density of state N{e) is given by the imaginary 
part of the real frequency Eilenberger function g(e) ob- 
tained by analytic continuation of the thermodynamic 
function g{uJn) from the imaginary Matsubara axis onto 
the the real energy axis We consider here the limit 
T = for which this procedure yieldsi^i^^: 



iy(e) = Im- 



2uvT^e 



tanx 



(40) 



where i/(e) — N{e)/N{0) is the normalized density of 
states, and x is the solution of the real-frequency version 
of Eq. (j2ip in which the substitution aj„ — —ie is made: 



A sin XT 

UVT — X 



uvT cot X ^ 1 



1 

4 



(41) 



Eqs. (gni) and (gl]) allow us to obtain N{e) using A{u) 
and Us calculated in the previous sections with the use 
of the thermodynamic Eilenberger equations. 

Figure [TUlfq) shows the evolution of i'{e) calculated for 
a clean superconductor (a — 0) as the condensate wave 
vector u increases. For u = 0, we recover the BCS density 
of states z^(e) = e/^/e^ — Ag for e > Aq and ^{e) — 
for e < Aq. In the presence of current, the singularity 
in j/(e) at e = Aq disappears, and the energy gap in 
the spectrum eg{a, u) defined as the maximum energy at 
which v{eg) vanishes, becomes smaller than Aq. Here the 
quasiparticle gap eg{u) decreases as u increases. At the 
critical value u = Us corresponding to the superheating 
field, the density of state at the Fermi surface in the 
gapless state equals j/(0) — 0.243. 

The effect of impurities on t^(e,M) for a moderately 
dirty superconductor with a ~ 3.6 is shown in Fig. [TU] 
(6). For M = 0, nonmagnetic impurities do not change 
the BCS density of states, in accordance with the Ander- 
son theorem. However, z^(e) in a current-carrying super- 
conductor with impurities begins to differ markedly from 
h'(e) in a clean superconductor. As it is evident from 
Fig. [TO] (a) and (5), impurities not only smear the cusps 
in h'(e) characteristic of the clean limit but also reduce 
A^(0) a,t u = Us, eventually restoring the gapped state at 
H = Hs where i/(e) = for e < eg{a). For the particular 
case shown in Fig. [TO] (6), our calculations give es(3.6) = 
0.211A(ms) = 0.169Aoo. For 1, the gap approaches 
the hmiting value = 0.410A(us) = 0.323Aoo.^-- 

The different behaviors of ^{e) in the clean and dirty 
limits shown in Fig. [TO] suggest that at the superheating 
field a quasiparticle gap es(a) appears as a superconduc- 
tor gets dirtier. This is illustrated by Fig. [TT] which shows 
the evolution of vie, Ug) at the depairing momentum as a 
increases. One can clearly see the transition from a gap- 
less to a gapped state induced by nonmagnetic impurity 
scattering, giving, for example, = O.ITAqo at a = 1. 
Therefore, at w = u^, a quasiparticle gap £^(0;) opens at 
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FIG. 10: Density of state v(e) in the current-carrying states 
for: (a) clean limit, a = 0; (b) moderate dirty limit, a = 3.6. 
Here A is understood as A(m) at T = for a given a. The red 
lines show v(e) at u = Us- Solid lines show i^(e) for u < Us, 
while dashed lines correspond to u > Us- For a — 0, the 
gap in the spectrum closed at u > 0.970us. For a = 3.6, we 
obtained es = 0.211A. 



a > ac, where the critical scattering rate ac calculated 
numerically from Eqs. (f40|) and (I4ip is: 

ac = 0.36 (42) 

The calculated quasiparticle gap (a) at m = is shown 
in Fig. 1121 Here es (a) monotonically increases as a in- 
creases above a > ac, approaching eg (00) « 0.323 Ago at 
a — ^ CX3. The dependence of eg on a can be approximated 
by the formula 

eg = 0.566Aoo(tan~^(0.626a + 1.345) - l) (43) 

to the accuracy better than 1.2%. 



V. NONLINEAR SCREENING AND Hs IN 
MULTILAYERS 

In this section we use the results obtained above to 
calculate screening of a magnetic field parallel to a mul- 
tilayer consisting of alternating superconducting (S) and 
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FIG. 11: Density of state v(e) at u = for different scatter- 
ing rates a. Here A = A(a) at T = and u = Us. 
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FIG. 12: Gap in tfie quasiparticle spectrum at tfie depairing 
momemtum u = Ms as a function of the nonmagnetic scat- 
tering rate a. The gap opens at a > 0.36. (Insert) The 
behavior of es{a) at small a. 



can be written in terms of the y-component u{x): 

j-^ = —J[u). 44 

We assume specular scattering of quasiparticles at the 
S-I interfaces, so there is no suppression of the order pa- 
rameter due to surface scattering in the S-layers. In this 
case J(u) is nearly uniform across each S-layer, a slight 
decreases of u{x) as x increases resulting from the Lon- 
don screening over the penetration length X ^ d. The 
solution of Eq. (|44l) gives the distribution of magnetic 
field B{x) = {(j)o/2TT)dxu{x) across the multilayer. 

The boundary conditions to Eq. (j44|) are as follows: 
dxu{0) — 2T:H/(j)o at a; = 0, u{oo) = at a; — )• oo, con- 
stant B in the Tlayers and the continuity of u{x) at the 
S-I interfaces. Shown in Fig. [13] is the distribution of 
H{x) calculated from Eq. (|44|) in which J{u) given by 
Eq. was calculated using the solution of the Eilen- 
berger equations for a = 1 and T = Q.bTc. As a compar- 
ison we also show H{x) (dashed line) calculated by solv- 
ing the linear London equation X^dxxU — u = for the 
same parameters. One can see that the Eilenberger the- 
ory which takes current pairbreaking into account, gives 
slower penetration of the magnetic field in the first few 
layers at the surface as compared to the London model 
which disregards current pairbreaking effects. 

1 



o 



0.8 
0.6 
0.4 
0.2 h 



insulating (I) layers, as shown in Fig. [T] Such multi- 
layers stabilize the Meissner state against penetration of 
vortices up to the superheating field of the material of 
the thin S-layers with c? ^ A for which the parallel lower 
critical magnetic field Hci = {2(j)o/Trd'^) ln(d/^) is greatly 
increased. In turn, the suppression of perpendicular 
currents in thin S-layers by non-conducting I-layers also 
suppresses the pairbreaking instability at the finite wave 
vectors k along the surface)^ which initiates penetration 
of vortices. As a result, the superheating field Hs in thin 
film multilayers is defined by the condition that J{u) in 
the first layer reaches the depairing current density. 

Penetration of magnetic field is described by the 
Maxwell equation V x V x A = AttJ / c, where the super- 
current J(u) depends on the gauge invariant phase gradi- 
ent u = Wip + 2TTA/(j)Q as shown in Fig. [2] (a). Here (p{r) 
is the phase of the order parameter. In the Meissner state 
B = V X A = ((/)o/27r)V x u, the Maxwell equation for 
the planar multilayer geometry in which u = (O, u{x), O) 



I \ \ \ \ \ 1 

0.5 1 1.5 2 2.5 3 

x/X 

FIG. 13: Distribution of the magnetic field Il{x) in a multi- 
layer obtained by numerical solution of Eq. (1441) for a — 1 
and T = O.STc (solid line). The result of the London model 
for the same parameters is shown by the dashed lines where 
A is the penetration depth for the material of the S layer. 

A curious manifestation of pairbreaking effects in the 
Meissner magnetization of moderately clean supercon- 
ductors occurs due to the transition from the gapped 
to the gapless state as the field H increases above Hg 
defined by Eq. Consider for example the magne- 

tization m{H) of a long cylinder of radius i? ^ A in a 
parallel magnetic field: 

1 

m = - / Jivs)dr (45) 
c Jo 
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where we use the dependence of J(us) on the superfluid 
velocity Vg = VcH{x)/Hg for T in the clean limit 
given by Eq. (|27|) . For H < Hg, the dc magnetization 
mo{H) = —H/Att corresponds to the linear Meissner ef- 
fect in the entire field region < H < Hg. However, 
as H exceeds Hg the magnetization exhibits singulari- 
ties in higher order derivatives indicating a phase transi- 
tion from the gapped to the gapless state. Indeed, let H 
be slightly above Hg so that m[H) = mo{H) + ma{H) 
where mo{H) — —H/Att is the ideal Meissner magneti- 
zation, and ma{H) is the nonlinear contribution due to 
the second term in Eq. (P7)) : 



c 



dx 



(46) 



where w{x) = H{x)/Hg, Hg = 47rJcA/c, H{x) H - 
Hx/\ is the field profile at the surface a; = 0, and L — 
[H — Hg)X/Hg <C A is the depth of the gapless layer. 
Equation (|46l) at e = {H — Hg)/Hg <C 1 then yields nia = 
(\/2/57r)iJg£^/2, which means a square-root singularity in 
the third derivative of ma{H) at H = Hg + 0: 



dm 



4ttV2H^ 



H 



-1/2 



(47) 



The singularity in the thermodynamic quantity 
d^m/dH^ implies the singular discontinuity in the 
forth derivative of the thermodynamic potential, in- 
dicating a field-induced forth order phase transition. 
Although Eq. (|47p was obtained for the clean limit 
a — 0, this singularity remains in the moderately clean 
limit as well if the scattering rate a < 0.36 is smaller 
than the critical value given by Eq. (|42|) so the transition 
from the gapped to the gapless state occurs below the 
superheating field H = Hg < Hg. Experimentally this 
rather weak singularity may be smeared out by local 
inhomogeneities of impurity concentration that result in 
a distribution of the local fields Hg{r). 



VI. DISCUSSION 

In this paper, we use the Eilenberger equations to cal- 
culate the effect of nonmagnetic and magnetic impuri- 
ties on the superheating field of type-II superconduc- 
tors. Unlike magnetic impurities which strongly suppress 
Hs, nonmagnetic impurities affect Hg weakly, although 
they can cause a nonmonotonic dependence of Hg on the 
scattering rate a at low temperatures. For instance, at 
T = 0, nonmagnetic impurities can increase Hs by ~ 4% 
at a ~ 0.6 as compared to the clean limit. As the scatter- 
ing rate a further increases, Hs{a) decreases and levels 
off at Hs{oo) sa Hs{Q), as shown in Fig. [S) This decrease 
of ffs at a ^ 1 is consistent with the decrease of the 
gap parameter As at the pairbreaking momentum u = Us 
shown in Fig.|31 Yet, unlike the decrease of the depairing 
current density Jc oc a"^/^ in the dirty limitj^S our re- 
sults show that the superheating field roughly scales like 



He even at low temperatures, so Hg is weakly affected by 
nonmagnetic impurity scattering. 

Our results obtained in the limit of k ^ 1 give the 
lower bound of Hs{T). The effect of finite n increases 
Hs{T) since the condition Vs{Q) — Vc is no longer suf- 
ficient to cause the instability of the Meissner state in 
the surface layer of thickness f where the superfluid ve- 
locities Vs{x) ~ (1 — x/X)vs{0) decreases below Vc be- 
cause of the London screening. The GL calculations^"— lii 
have shown that the finite-K effects increase Hs{k) ~ 
(l -I- 0.7K~^/2)i?s(oo), giving a correction ~ 7 — 16% as 
compared to Hs{oo) at k — ^ oo for k = 20 — 100. These 
effects also make Hs(k,) dependent on a since K(a) in- 
creases as a increases approaching k ^ ano in the dirty 
limit. Thus, Hs{T) generally decreases as the scatter- 
ing rate a increases although this effect is comparatively 
weak for high-K materials. Addressing the dependence of 
Hg on a due to the finite-K effects at low temperatures 
requires the instability analysis of the Meissner state with 
respect to 2D perturbations of Sf{x,y) and SJ{x,y) de- 
scribed by the linearized Eilenberger equations. 

The nonmonotonic dependence of H^ on a obtained in 
this work results from interplay of Meissner currents and 
impurity scattering, and their effect on the quasiparti- 
cle density of states A^(e, J). Our calculations revealed 
the disorder-induced transition from the gapless to the 
gapped state at H = Hg, which can have important im- 
plications for the low-temperature surface resistance Rg 
at high RF fields H{t) ^ He- The BCS surface resistance 
Rs at small RF fields H <^ HcT/Aq, low frequencies 
(w < Ao) and T < is given by2ii25 



Rs =uj — - — exp 



Ao 
T 



(48) 



where the factor A depends on the mean free path £ and 
(weakly) on the RF frequency to. The main Boltzmann 
factor exp(— Aq/T) accounts for the exponentially small 
density of thermally-activated quasiparticles due to the 
zero density of states N{e) for the energies e < Aq. In 
the presence of Meissner current, the quasiparticle gap 
£g{J) shifts to smaller energies, giving rise to a highly 
nonlinear dependence of the surface resistance on the RF 
field amplitude at low temperatures: 



Rs oc exp 



AH) 
T 



(49) 



As the field increases, Rs{T,H) becomes essentially de- 
pendent on H if the field-induced change of eg{H) is of 
the order of In the clean limit for which eg{H) — 
(1 — H/Hg)Ao, this condition takes the form 



H > TH,/Ao. 



(50) 



Therefore, the dc superheating field Hs in the clean limit 
has no direct relevance to the maximum RF magnetic 
field at which the Meissner state can exist. Moreover, 
even at fields H^ smaller than Hg < He, the RF field 
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starts generating quasiparticles as the gap eg{H) becomes 
smaller than the RF frequency: eg{H^) — hoj. Hence 

Hg-H^> hioHg/Ao (51) 

As the magnitude of H{t) approaches H^^, the surface 
resistance increases strongly, becoming of the order of 
Rs{T) in the normal state. Calculation of the non- 
linear surface impedance requires solving equations for 
the nonequilibrium Keldysh functions that take into ac- 
count not only current pairbreaking, but also the effect 
of RF field on the quasiparticle distribution function de- 
termined by collisions of electrons with impurities and 
phonons^. We will not discuss here this complex prob- 
lem in detail, but only make a few qualitative remarks 
based on our solutions of the dc Eilenberger equations 
that capture the essential effect of impurities on the 
quasiparticle density of states at high RF fields. 

Nonmagnetic impurities can strongly affect the field 
dependence of Rg because they reduce by only a few 
percentages but restore the gap eg in the quasiparticle 
spectrum a.t H = Hs, as shown in Figs [TT] and [T^l In the 
dirty limit, a > 1, the gap eg ~ O.SAqo may therefore be 



big enough to ensure both the exponentially small surface 
resistance in Eq. at low temperatures and the lack of 
quasiparticles generated by the RF field with uj < eg. As 
a result, nonmagnetic impurities can drastically reduce 
the field-induced increase of Rs as compared to the clean 
limit. Moreover, the dc superheating field Hg « 0.83i?c 
at a > 1 can now be regarded as a true maximum field 
amplitude at which the Meissner state can survive under 
low frequency RF fields. This conclusion may be essen- 
tial for the materials optimization for superconducting 
cavities used in particle accelerators. 
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